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Abstract 



Fluctuations in the cosmic microwave background (CMB) temperature are 
being studied with ever increasing precision. Two competing types of theo- 
ries might describe the origins of these fluctuations: "inflation" and "defects". 
Here we show how the differences between these two scenarios can give rise 
to striking signatures in the microwave fluctuations on small scales, assum- 
ing a standard recombination history. These should enable high resolution 
measurements of CMB anisotropies to distinguish between these two broad 
classes of theories, independent of the precise details of each. 

PACS Numbers : 98.80.Cq, 95.35+d 
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Modern experiments are producing a growing body of data on the fluctuations in the 
cosmic microwave background (CMB) Q. The origin of these fluctuations may be due to 
defects |fJQ] such as cosmic strings or textures, or an early period of cosmic inflation M. Our 
understanding of both pictures is sufficiently incomplete to allow a wide range of predictions 
from each []5]||. Recent calculations of degree scale anisotropies from defects have focused 
on the case of a reionized thermal history J7|. Here, we consider a standard recombination 
universe, and show how fundamental differences between the defect and inflationary scenarios 
can leave sufficiently different signals as to allow the high resolution CMB measurements to 
distinguish between them. 

The history of radiation can be separated into three epochs: at very early times (the 
"tight coupling" epoch) photons were tightly coupled to baryonic matter and together they 
behaved like a relativistic fluid in which pressure waves propagated at c/y/3. As the uni- 
verse expanded this coupling weakened, eventually producing today's "free-streaming" epoch 
where photon-matter interactions are negligible. These two epochs are linked by a "damp- 
ing epoch" during which dissipation can diffuse perturbations [PHT0|]. We will focus our 
attention on the first of these epochs. 

Due to their quantum nature, inflationary models predict an ensemble of possible states 
for the universe, and the ensemble of sub-horizon pressure waves of the photon-baryon 
fluid in an inflationary model are depicted in Figure 1 (for fixed wavelength). Note that the 
ensemble is "phase focused" , so that the entire ensemble achieves zero amplitude at the same 
instant in time |TT|. These moments of zero amplitude occur at different times for different 



wavenumbers, and the result is the striking oscillatory behavior in the power spectrum at 
the end of the tight-coupling epoch (Figure 2, solid curve). 

In contrast, cosmic defects will drive sub-horizon pressure waves in a random manner 
which tends to destroy the phase focusing (see Figure 3). In many realistic cases (such as 
cosmic strings) we expect this effect to prevent an oscillatory power spectrum from emerging 
(Figure 2, dashed curve). Figure 4 illustrates how these differences translate into observable 
features in the microwave sky. For other cases (such as cosmic texture) the phase decoherence 
of the evolving defect network is less effective, but still the oscillations are of sufficiently 



different origin to leave a telltale signature [12,13 



We now treat each of these points in more technical detail. We evolve the radiation 
and matter fluids from some early epoch up to the time of last scattering and examine 
the processes which create or destroy the characteristic oscillations at decoupling. The 
calculation of the angular power spectrum (Figure 4) involves additional physics and is 
treated in a companion letter [13[ . 



We work in the synchronous gauge, and use the variables of reference |15[| , extended to 
include a baryonic component. In Fourier space these are 

r o = e D + -^-(-)Vs[i + J R] (i) 

5 C = A'kG- (r 00 - OO ) - - ( ^ c + 2 [1 + R] Q r ) 5 C 

- -2 [1 + R] Q r s (2) 



a 
R 



1 + R 



s - c 2 s k 2 (s + 5 C ) (3) 
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FIG. 1. Perturbations from Inflation: Evolution of two different modes during the tight coupling 
era. While in (a) elements of the ensemble have non-zero values at r]*, in (b), all members of the 
ensemble will go to zero at the final time (77*), due to the fixed phase of oscillation set by the 
"growing solution" initial conditions. The y-axis is in arbitrary units. 

Here is the pseudo-stress tensor, Q D = di<d 0i , is the defect stress energy, a is the 
cosmic scale factor, G is Newton's constant, 5x is the density contrast and fix is the mean 
energy density over critical density of species X (X = r for relativistic matter, c for cold 
matter, B for baryonic matter), s = |<5 r — 5 C ,R= \pb/ Pr, Pb and p r are the mean densities 
in baryonic and relativistic matter respectively, c s is the speed of sound and k is the comoving 
wavenumber. The dot denotes the conformal time derivative d v . 

The evolution described by these equations is very different for the two classes of theories 
we are considering. Consider first perturbations of a primordial origin, such as inflation, 
with no defects (0 MJ , = 0). At re-heating, all modes of interest are outside the horizon, 
with random phases and amplitudes chosen from some Gaussian distribution. On these 
super-horizon scales the radiation overdensity has one growing and one decaying solution 
(for example, in an inflationary model these behave as rf and r/~ 2 in the radiation-dominated 
epoch). Just before horizon crossing these perturbations from inflation (and, more generally, 
most perturbations set by "initial conditions") will be dominated by the growing solution. 
If horizon crossing occurs in the tight-coupling epoch the growing solution will match onto 
a temporally oscillating solution, a pressure wave in the relativistic fluid. 

The matching of a super-horizon pure growing mode to the sub-horizon pressure wave 
forces the wave to have a unique temporal phase (for a given mode, k) across the entire 
ensemble. For example, in Figure la we see the evolution of several members of the ensemble 
for a specific k. For each realization (or member of the ensemble) the pressure wave has 
exactly the same phase. In Figure lb we pick another value of k for which, regardless of 
which element of the ensemble we take, all S r (k) lie on a node at 77*, the end of the tight- 
coupling epoch. Figure 2 shows the resulting r.m.s. of 5 r (k,r]+) (solid curve) for a range of 
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FIG. 2. The r.m.s. value of 5 r evaluated at decoupling (?7*) for inflation (solid) and cosmic 
strings (dashed). For this figure we use J%) = (1 + 2(kr]) 2 ))~ 1 , To = 1/(1 + (2vr/A;7?) 2 ) 2 , and 
rj c = 77/(1 + krj). 



k. The super- horizon growing condition gives $,.(77*) = only for particular values of k, and 
the net result is the oscillatory spectrum above. 

For cosmic defects the story is very different. In this case equations ([j])-(JD are sourced 
by the functions Qoo(k,r]) and Qoik,^) (see for example [0). Defects evolve in a highly 
non-linear manner, and nothing short of a large, three dimensional numerical simulation 



can be expected to give a M ^(/c, 77) correct in every detail |I1||l"6|j . Here, we resort to a 



simplification which we expect to be a good approximation to the true defect evolution, and 
which certainly allows us to illustrate our main points. 

We are interested in evaluating the radiation power spectrum at decoupling and therefore 
need not concern ourselves with detailed knowledge of the statistical properties of the defect 
network. It suffices to model correctly the one- and two-point correlation functions of the 
defect stress-energy. In Fourier space this means modelling the equal and unequal time 
correlation functions of Goo and 0/> (We will assume zero cross-correlation between the two 
source functions, a freedom we can take without violating defect stress-energy conservation.) 
Let us construct functional forms for the equal time correlations first: Scaling of the network 
will restrict (|Ooo(&, v)\ 2 ) = J~oo(kr]), where causality ensures that JF 00 oc (kr])° on large 
scales (kr] <C 1). The covariant conservation of the defect energy- momentum tensor forces 
(|@r>(&, v)\ 2 ) = ^^D{krj), where JF D oc {krj) A on large scales. Depending on the type of 
defect, Too and Tr> will have different behavior on small scales; for cosmic strings, JF 00 oc 
(kr])~ 2 on scales smaller than a fraction of the horizon [17], while for textures the turn-over 



occurs on scales of order the horizon, and has a much steeper (kf])~ 4 behavior ||18|| . (The 
caption to Fig 2 gives one choice for these functions which models cosmic strings.) 

Next we note that the urTequal time correlation functions (such as (Qoo(k, r])Q 00 (k, r] + 
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FIG. 3. Perturbations from defects: Evolution of 6 r (k) and the corresponding source 0oo during 
the tight coupling era (0£> is not shown). Two members of the ensemble are shown, with matching 
line types. Due to the randomness of the source, the ensemble includes solutions with a wide range 
of values at 77*. Unlike the inflationary case (Figure 1) the phase of the temporal oscillations is not 
fixed. The y axis is in arbitrary units, and the source models are the same as for figure 2. The 
factor rja/a allows one to judge the relative importance (over time) of the 0oo term in Eqn 2. 

A77))) will go to zero at some finite value of A77 characterised by a "coherence time" , i] c (k, if). 
The lack of long time coherence is due to the fact that each mode is coupled in a highly 
non-linear way to all the others. To a first approximation this coupling may be thought of 
as producing "random" kicks to a given mode, driving the unequal time correlation function 
to zero. 

Using only the equal time correlation function and the coherence time we construct an 
ensemble of realizations of ©oo^ 7 ?) (and similarly Ojj(A;, 77)) in the following way: At the 
initial time 77, we choose a value for Qoo(k,r)i) from a Gaussian distribution with variance 
given by (|Ooo(&, Vi)\ 2 )- We next choose a time step from a distribution with mean value r/ c , 
and choose a value for Boo a t the new time, again from a distribution with variance given 
by this equal time correlation function at the new time. We continue this process until the 
final time is passed, and construct a smooth function using a cubic spline. The process is 
repeated many times to produce the ensemble of source histories. 

All that remains is to choose initial conditions for the remaining variables of equations 
(ID - ®) an d here causality comes in. The formation of defects is assumed to occur against 
a completely homogeneous and isotropic background. It occurs in a finite time, and so 
is limited by causality to only move matter a finite distance (of order the horizon at the 
defect-producing phase transition). This translates into strong k 2 suppression of Too on 
large scales (see p!5|J19|j20| ) . To reflect this we choose initial conditions with r o(?7i) = for 
all k. For a given initial defect source, Too = implies a "compensating" perturbation in 
the radiation and dark matter on large scales. The advantage of this formulation of the 
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FIG. 4. Angular power spectrum of temperature fluctuations generated by cosmic strings 
(dashed) and arising from a typical model of scale invariant primordial fluctuations (solid) in arbi- 
trary units. The all-sky temperature maps are decomposed into spherical harmonics 4? = a^^m 
from which one defines the angular power spectrum as C\ = 2i^rZ)-«l a ml 2 - The shape of the string 
curve for / ~ 100 (and thus the height of the peak) is very sensitive to existing uncertainties in 
string networks. We show only the scalar contribution, in arbitrary units. The string curve is from 
[13] where we use an extended Hu-Sugiyama formalism. 



perturbation equations is that with such a choice of initial conditions, the compensating 
physical constraints are built in and require no further attention. 

Figure 3 shows the evolution of the 5 r in a cosmic string scenario, with the corresponding 
source evolution for each realization. The model is sufficiently incoherent to produce random 
phases in the pressure waves. Cosmic string networks are not well enough understood to 
dictate precisely the functions J-d,oo, an d Vc(v)- We have explored a range of choices for 
these functions suggested by string simulations, and have found that they all have enough 
small scale power (and a small enough typical coherence time, to be well inside the 
decoherence-dominated regime pip . Thus in all these cases the power spectrum of the 



radiation at decoupling (the dashed curve in Figure 2) does not show the oscillatory features 
found in the inflationary curve. 

A particular high-coherence limit (taking = ((|@(/c, rj) | 2 )) 1//2 ) has been explored by 
others in the context of the cosmic textures [|T2|, 1S|] - The models they used to calculate the 



microwave sky include none of the decohering effects we discuss here and, naturally enough, 



oscillations are produced. Crittenden and Turok |12| motivate the high- coherence limit by 
showing various correlation functions (including ones similar to Figure 2) measured in large 



1 Even the rather conservative choice r\ c = rj does not give oscillations. 
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numerical simulations. The texture simulations show significant oscillations, indicating that 
a high level of coherence is present. 

After the damping era, the CMB photons stream towards us and inhomogeneities in 
S r on different scales k will be converted into anisotropics AT/T seen today on different 
angular scales. Gravitational perturbations between the surface of last scattering and now 
will add on further anisotropics, but only on angular scales larger than the ones affected 
by the pre-?]* physics. Hence the I > 100 angular power spectrum will project into our sky 
the 5 r (k,r]+) peak structure, the so-called "Doppler peaks" (softened somewhat by the effect 
of the radiation fluid's velocity). In [T3| we account for all these effects by extending the 
formalism of Hu and Sugiyama |22| to include defects. 

As an illustration, Figure 4 gives the angular power spectrum of the brightness of the 
microwave sky for the standard CDM inflationary model and one of the string-based models 
considered in ||13|| . The string curve exhibits complete suppression of the secondary Doppler 
peaks, a signal which should be easily resolved by high resolution experiments (such as 
COBRAS/SAMBA p3|). Although the degree of secondary peak suppression can depend 



on the details of the defect model ||13|| , we have found that realistic string models all appear 
to achieve complete suppression We emphasize that the shape of the string curve at 
low I and the peak position are still subject to uncertainty [|T3|,^T|]. This is particularly true 
of the peak height, which is uncertain within an order of magnitude. 

We have seen how the high degree of coherence which is present in inflationary scenarios 
(and some defect scenarios) leads to phase focusing of the sub-horizon pressure waves. This 
effect leads to "secondary Doppler peaks" in the angular power spectrum. For defect scenar- 
ios the random effects of non-linear defect evolution tend to decohere this focusing. Models 
representing cosmic strings exhibit a striking suppression of the secondary Doppler peaks in 
the angular power spectrum. Even in cases which are not so extreme, we expect that the 
decohering effects will lead to distinct observable signals in the angular power spectrum. 
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